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Abstract —Due to the increasing concern for greenhonse gas 
emissions and fossil fnel security, electric vehicles (EVs) have 
attracted much attention in recent years. EVs can aggregate 
together constituting the vehicle-to-grid system. Coordination of 
EVs is beneficial to the power system in many ways. In this 
paper, we formulate a novel large-scale EV charging problem 
with energy trading in order to maximize the aggregator profit. 
This problem is non-convex and can be solved with a centralized 
iterative approach. To overcome the computation complexity 
issue brought by the non-convexity, we develop a distributed 
optimization-based heuristic. To evaluate our proposed approach, 
a modified IEEE 118 bus testing system is employed with 10 
aggregators serving 30 000 EVs. The simulation results Indicate 
that our proposed distributed heuristic with energy trading can 
effectively increase the total profit of aggregators. In addition, 
the proposed distributed optimization-based heuristic strategy 
can achieve near-optimal performance. 

I. Introduction 

With the increasing concern for greenhouse gas emissions, 
electric vehicles (EVs) are expected to reach a signihcant 
market share in the near future. With the emerging vehicle- 
to-grid (V2G) technologies, EVs can potentially help alleviate 
the security concern for the supply of fossil fuels and mitigate 
the power network instability caused by the highly penetrated 
intermittent renewable energy generations Q. However, the 
network security and economic operation can be signihcantly 
adversely influenced by the uncoordinated charging behaviors 
of a large number of EVs Q- By efficiently utilizing the 
system capacity, coordinated charging strategies can reduce the 
possible adverse impacts on the power system |j^. In addition, 
other merits can also be obtained, such as reducing the total 
operational cost and mitigating the variability of the renewable 
energy sources Q. Hence, coordination of EVs is beneficial 
to the power system in many different perspectives. 

A large population of EVs can be clustered into aggre¬ 
gators to facilitate the coordination of charging behaviors. 
Coordinated EV charging within one aggregator has been 
studied extensively in recent years. Previous research generally 
focused on the grid structure with aggregators managing a 
large number of EVs. Based on where the charging decisions 
are made, the methodologies can be classified into two types: 
centralized and distributed methods. In the centralized meth¬ 
ods, EVs need to send their parking and battery information 
to the associated aggregator, which then decides how each 
individual EV should be charged in a centralized manner. 


The decision is mostly driven by maximizing the aggregator 
profit minimizing EV owners’ power cost 0, 0, and 

achieving power balance 0. On the other hand, aggregators in 
the distributed methods send power pricing scheme and related 
information to EVs. The EVs then utilize their own knowledge 
to establish their charging plans, which are later delivered to 
the aggregator. Usually this process repeats until an agreement 
or equilibrium is reached. As the charging strategies are 
developed by the EV, most existing efforts employ the energy 
costs incurred for the EV owners as the performance metric 
0- User utilization ratio Gg and power balance im are also 
studied as metrics. In general, all these mainly focus on the 
coordinated charging behavior within a single aggregator. 

There is few efforts studying the interactions among mul¬ 
tiple EV aggregators, e.g., 0, nig. One common feature 
is the utilization of a hierarchical architecture where ag¬ 
gregators coordinate with others under the control of the 
system operator. While this architecture efficiently reduces 
the computational burden of developing coordinated charging 
plans, direct interactions among aggregators are generally not 
considered. 

In this paper, we focuses on maximizing the aggregators’ 
total profit through coordinated EV charging when considering 
multiple aggregators connected to the grid. Different from 
the previous work, we investigate the possibility of achieving 
Inter-aggregator Energy Trading (lET) to further exploit profits 
from trading energy directly among aggregators. There exists 
previous research investigating the possible energy trades 
among EVs, e.g., GD and Gg, and they demonstrated im¬ 
proved total costs of EVs when compared with those models 
without energy trading. Therefore it is likely that energy 
trading among aggregators is more profitable with direct 
aggregator interactions. In addition, we consider the Locational 
Marginal Pricing (LMP) in this work to manage the system 
congestion cost to further increase the profit. 

The remainder of this paper is organized as follows. Section 
II presents the system framework of our proposed approach. 
We formulate the aggregator profit maximization problem and 
propose a centralized coordinated charging strategy in Section 
III. Section IV demonstrates an optimization-based heuristic 
approach for the problem. In Section V we provide a case 
study to illustrate the performance of our developed approach. 
Einally we conclude this paper in Section VI. 


II. System Architecture 

In this work, we adopt a typical hierarchical EV charging 
architecture which is composed of three components: 

system operator (SO), aggregators, and plug-in EVs. The 
main purpose of introducing aggregators is to simplify the 
system model and unbundle the electric loads, e.g., EVs, 
from the power network infrastructure HD- By collecting 
the charging information from the clustered EVs and pricing 
information from SO, each aggregator can develop its own 
optimal charging schedule such that the aggregators’ profits 
are maximized while the various EV charging requirements 
are all satisfied. SO manipulates the real-time power price to 
balance power consumption of each bus in the power grid to 
alleviate power congestion. 

We consider EVs as dispatchable power appliances whose 
charging rates can be adjusted by the corresponding aggrega¬ 
tors by taking the system requirements and economic consid¬ 
eration into account |Tg. Similar to GZl’ upon the arrival of 
an EV, we assume that its battery capacity and current state of 
charge (SoC) can be obtained by its aggregator via appropriate 
vehicular communication techniques. In addition, the departure 
time of the EV is also assumed to be available to its aggregator. 
We allow EVs to perform early departure if needed. In such 
case, a penalty can be imposed to the EV owner p8) . 

In our proposed framework, EVs are categorized into two 
classes based on their respective capabilities: uni-directional, 
i.e., charging from the grid to vehicles, or bi-directional, i.e., 
both charging from and discharging to the grid. Enabling 
power flow from EVs to the grid can effectively increase the 
profit by exploiting the power price fluctuations 0,0, ng. 
To encourage V2G bi-direction operation, aggregators may 
compensate the EV owners’ battery aging loss by offering 
lower charging cost or free battery replacement plan | [20| . 

The objective of lET is to effectively utilize the discharging 
power from aggregators for the benefit of aggregator profit. 
The aggregators typically purchase power from SO for charg¬ 
ing at a price, and sell power to SO at a lower price. Mean¬ 
while, instead of trading with SO, the discharging aggregators 
can trade with other charging aggregators at a better price, 
between the SO power buying and selling prices. Eor instance, 
suppose that Aggregator A1 is charging power Pi > 0 from 
the grid with power purchasing price and Aggregator A2 is 
discharging power P 2 < 0 to the grid with power selling price 
Then the total cost for the aggregators is PiCf'+P 2 C 2 ^. 
Suppose that A2 desires to sell > 0 to A1 at The 

power cost for A1 becomes _ptrade^( 2 ^ch_|_ptrade(^trade 
A2’s is (P 2 -f pf^de)(^dch _ ptradej^trade^ rendering a total cost 
at PiCf' + P 2 Cf^ + Therefore the trades 

between A1 and A2 can decrease the cost, thus increase the 
aggregators’ total profit when C'^'^ > Cf'. 

III. Aggregator Profit Maximization Problem 

In this section, we formulate the aggregator profit max¬ 
imization problem with lET We employ model predictive 
control to develop charging control strategies of the immediate 
time slot, taking future EV departures and power price profile 


into account. Specifically, the problem is defined over a finite 
time-horizon T = {tq\tq = fo = 0, lu'' 

Information in the current time slot to is considered accurate 
while the future information may be imprecise. Since the 
solution of the problem, i.e., the overall charging strategy, 
should be jointly optimal over the entire T, the implemented 
charging strategy of the current time slot will also contribute 
to the total profit maximization over T, instead of one single 
time slot. The problem shall be solved whenever EVs require 
charging operations, rendering the process online. All symbols 
used in this paper and their meanings are listed in Table |I] 

In the problem, a set of m aggregators 

A = {Ai, A 2 , - ■ ■ ,Ai,--- ,Am} are considered. 
Aggregator Ai serves EVs at a time denoted by 

Vi = Vj, 2 , • • • , Vj,riJ. The objective function 

for time slot t is mathematically formulated as follows^ 

E + E E 

AiGAVijGVi Ai^A VijGVi 
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Chaiging Income Penalty Income 

- E ( E (1) 

Ai^A Vij^Vi 


Marginal Energy Cost 


where 





if Ey..,6V. 

otherwise, 


(2) 


which means that the charging price is employed if the aggre¬ 
gator draws power from the grid. Otherwise the discharging 
price is used. The income of aggregators is composed of two 
parts: charging and penalty incomes. The first term in Eq. 0 
is the charging fee inmosed on the EVs, which is arbitrarily 
formulated as followso 


Ci 


o'" — AC'^ min{|7ijj/r'^, 1} bidirectional charging, 

— AC'^ min{|7ijj/r^, 1} unidirectional charging, 

(3) 


where the bidirectional charging allows both charging and dis¬ 
charging behaviors, and the unidirectional charging supports 
charging only. The second term in Eq. 0 is the penalty in¬ 
come (parking fee) imposed on the EVs for any late departure. 
This income can be neglected in residential area scenarios. 
The third term in Eq. 0 is the cost of the consumed power 
purchased from SO. The cost is firstly generated using day- 
ahead energy prices, then taking LMP into consideration in 
optimization of the later time slots. If the charging power 
is negative, the corresponding aggregator will perform V2G 
energy selling operation and this term will become negative. 

Besides real time profit generation, aggregators can also 
utilize the dispatchable characteristics of EVs to delay the 
charging process to further increase the profit. This results 
in a mult-time-slot joint optimization as follows: 


*For the sake of simplicity, the symbol t for time is omitted in equations 
(l)-(3), (6), and (8)-(ll). when no confusion may be caused. 

^More complicated contract formulations are available in practice and can 
be easily incorporated into the proposed optimization problem. 






TABLE I 

Sets, Parameters, and Variables Used 


Parameter 

Description 

Pai'ameter 

Description. 

A 

Set of aggregators in the system. 

plme 

Maximum power flow on line 1. 

Ai 

2 -th aggregator in the system. 

/^mm /^max 

’^S 

Minimum and maximum generation output for the 

W 

Set of Electric Vehicles (EVs) served by Ai. 


generator on bus 5. 

Vi,3 

j-th EV served by Ai. 

t 

Current time slot. 

Vi,u 

u-th unidirectional charging EV served by Ai. 

Rt 

The total aggregator profit at t. 

T 

Time horizon of optimization. 

Pi 

Total charging power of Ai. 


Parking time of Vij. 

Pi,3 

Charging power of Vij. 


The current and the q-th future time slot. 

^i,j 

Binary indicator for whether Vij is still parked after 

gHiax 

The total number of future time slots. 


the registered departure time. 

At 

Duration of each time slot. 

^i,3 

SoC of Vij. 

■ dpi 

The registered departure time of Vij- 

^vmin 

Minimum reserved State-of-Chai‘ge (SoC) of Vij- 

Pt 

Power charging price from the power grid at Ai. 

Gs 

Generation output for the generator on bus s. 

^ach 

Power discharging price from the power grid at Ai. 

Ps 

Aggregator injection on bus s. 


The base contract price for bi-directional and uni¬ 
directional vehicle-to-grid (V2G) charging enabled EVs. 

Ci,j 

Charging fee of Vij when the registered parking time 
is 7i,j . 

AC'',ACC 

The marginal price decrease for bi-directional and 

^trade 

Power trading price among aggregators. 


uni-directional V2G charging enabled EVs. 

Ci 

Power trading price from the power grid at Ai. 

rpW rpU 

Parking time for minimal charging fee of bi-directional 

Dtrade 

Trade power of Ai. 


and uni-directional V2G charging enabled EVs. 

pspl 

Total available aggregator power supply in the trade 

r>ch odch 

Maximum chai'ging and discharging rate of Vij. 


at Ci. 

„ch „dch 
''i,j ’'iij 

Charging and discharging efficiency of Vij. 

Ddmd 
^ i 

Total available aggregator power demand in the trade 


Battery capacity of Vij. 

pcap 

at Ci. 

cmiri cmax 

Minimum and maximum State-of-Charge (SoC) of Vij. 

Trade capacity at Ci. 

5.1 

Required SoC on departure of Vij. 

Cs{Gs) 

Power generation cost function for the generator on bus 

N 

Number of buses in the power system. 


s when generating Gs power. 

L 

Number of lines in the power system. 

\tii 

Lagrangian multipliers. 

Ds 

Fi-s 

Inelastic load on bus s. 

Generation shift factor from bus s to line 1. 

_min _max 

Is ’ Is 

Lagrangian multipliers. 


maximize > Rt subject to 

p pirade ' ^ 

ter 

0 < Pi,u,t < Ptu,^t G r, 


- t - l)Jt; > (50 - S^,j,t)E,,,/nt 


Vi G r, 


p.,j,t/^tntj < (50 - G r, 

E 


= o,Vi G r, 






E 


i,t 


> o,Vi G r. 


(4a) 

(4b) 

(4c) 

(4d) 

(4e) 

(4f) 

(4g) 

(4h) 


Eq. ( |4b] i imposes rigid upper bounds for charging power of all 
EVs online. Eq. ( |4c| further limits the uni-directional charging 
EVs to perform charging operations only during the whole 
parking process. Eq. ( |4d| ) ensures that the current charging 
power of an EV is feasible only if the battery can be charged to 
the required SoC when the charging powers of all succeeding 
time slots are set to the maximum rate. This constraint is con¬ 
sidered when the optimization horizon T does not contain the 


departure time i.e., max{T} < Otherwise, 


an alternative constraint is considered instead of Eq. (4di, as; 




t&r 


(5) 


This constraint guarantees that EVs will be charged to their 
required SoC on departure. Other practical strategies to prevent 
insufficient charging, e.g., penalties on aggregators, can be 
easily implemented in our model. Eq. ( |4^ prevents the battery 
from being over charged by limiting the current charging 
rate. Eqs. ( [4d| and ( [4^ cooperate to manipulate the charging 
rates of EVs to satisfy the SoC constraints. Eq. ensures 
that the amount of energy bought and sold in the aggregator 
trading market are equal. Non-convex constraint Eq. ( |4^ 
prevents the aggregators, that are consuming power from the 
power grid, from trading energy to other aggregators, and 
vice versa. As and J^Vi gv, are independent and 

can be both positive and negative. Eq. imposes that 
the aggregators cannot sell more energy than they generate 
through discharging behaviors, and cannot buy more energy 
than they request to fulfill the EV charging demand. 

When given constant power prices Q’s, the optimization 
problem formulated in Eq. 0 can develop optimal EV 
charging profiles for maximizing aggregator profit. However, 
the optimized charging powers of each aggregator Pi = 
Thv GV change the power flow of the grid, 

resulting in changing C^’s. Therefore, we also employ an 
Optimal Power Elow (OPE) based LMP optimization together 
with Eq. 0 to develop C^’s subject to aggregators’ changing 
Pi values. 

Upon the receipt of charging power requests from the 
aggregators, SO performs OPE to route the power flow. In 
our proposed framework, an LMP strategy is employed for 











Fig. 1. Flow chart for the proposed centralized coordinated charging strategy 
for aggregator profit maximization. 



Fig. 2. Flow chart for the proposed optimization-based heuristic coordinated 
charging strategy for aggregator profit maximization. 

congestion control HD- As we focus on the real time power 
dispatch of the aggregators for EV charging, it is assumed 
that the Unit Commitment problem has been solved and 
all generators considered are online Q. The SO level OPE 
problem is formulated as a classical OPE problem | [22) , where 
the LMP of buses are optimized with the partial derivative of 
the Lagrangian of OPE (see for elaborations). 

Optimization problem Q and OPE-based LMP optimization 
together compose the aggregator proht maximization problem. 
This problem aims to develop optimal EV charging plans 
to maximize the aggregator proht. It can be solved in an 
iterative manner as depicted in Eig. Problem and OPE 
are optimally solved alternately until an equilibrium is reached, 
where a stable lET-enabled EV charging prohle is developed. 
At the beginning of each time slot, Eq. Q is solved to generate 
optimal charging strategies for all EVs in the system, using 
the day-ahead energy price. After the scheduled aggregated 
charging powers are optimized, they are reported to the SO 
for calculating the new LMP. Utilizing these new prices, the 
central controller solves Eq. 0 again, and this process repeats 
until an optimal charging strategy is developed. 

IV. Distributed Aggregator Heuristic 

The aggregator proht maximization problem develops op¬ 
timal EV charging prohles for maximizing aggregator proht. 


However, the problem can become intractable when there is a 
large number of EVs to schedule. As Eq. 0 is a non-convex 
problem, existing solvers may encounter efficiency issues 
when solving large problem instances. In order to address this 
problem, we propose a distributed heuristic algorithm to solve 
Eq. Q. The heuristic-enabled aggregator proht maximization 
strategy is shown in Eig. We basically divide the shaded 
step in Eig. [T] (which constitutes the non-convexity) into the 
three shaded sub-steps in Eig. They respectively correspond 
to the three steps of the proposed heuristic: aggregator proht 
optimization (shaded sub-step on the left in Eig. |^, power 
trading (top-right), and supply-demand balancing (middle- 
right). 

A. Aggregator Profit Optimization 

At the beginning of each optimization iteration, all online 
EVs send their vehicle and battery information, including the 
maximum charging rate, scheduled departure time, and bat¬ 
tery size and current SoC, to their corresponding aggregator. 
With this information, each aggregator maximizes its proht 
independently using a modihed formulation of Eq. 0: 

Ri,t = -f (1 — Cei,j)PijCiJ 

Vi.,ev, VijeVi 

- ( y]] ai,jP^,j - (6) 

VijeVi 

The major difference between Eqs. 0 and 0 lies in the 
introduction of the trading cost term Pf'“*‘^C'trade- is set 

to zero at the beginning of the iteration, and set to either a 
constant value or to gv based on the supply 

and demand equilibrium, which will be elaborated in Section 
IV-B. Therefore, Eq. Q is transformed into 

max^^ Ri t subject to ( |4bl l, ( |4c] l, ( |4dl l or 0, ®, (7) 
ter 

to maximize each aggregator’s proht. Eq. 0 considers all EV 
constraints, and energy trading constraints Eqs. @-(l^ are 
handled in the power trading heuristic algorithm. 

B. Power Trading Heuristic 

When all aggregators have hnished their proht maximization 
processes, they broadcast their Pi values to the others in 
the aggregator trading market in the form of bids. In this 
step, the aggregators utilize all the power bids to perform 
a modihed second-price auction | |23| . When submitting bids, 
each aggregator calculates its Pi from the result of Eq. 0^ 
and places a bid in the form of {Pi,Ci) pair. A positive Pi 
makes Ci = Cf', and a negative Pi makes Ci = The 

bids are then broadcast to the others and the auction starts in 
the aggregator trading market when all bids have been placed 
and announced. 

After all bids are generated, the auction evaluates the 
available charging (demand) and discharging (supply) power 
for trade at all possible trading prices whose values are 

selected from all Ci values in the bids: 

A^eAf AkdAfi'^ 


( 8 ) 




































where 

= {Afe e A\Pk <0,Ck< (9) 

4"“^ = {Ak G A\Pk >0,Ck> (10) 

The trading capacity for is accordingly calculated by 

PfP = min{-P’*P\ Pf ^^} X (11) 

Consequently, the value of is set to the optimal Ci value 
which yields the largest 

After determining the trading price, another round of opti¬ 
mization Eq. (j7]i is conducted to utilize the additional profit 
contributed by the trading behavior. The generated charging 
profile is considered final and reported to SO for LMP 
calculation. 

V. Case Studies 

We employ the IEEE 118 Bus system p4) to assess the 
profit maximization performance of our proposed approach. 
10 aggregators are installed in the system on Buses 7, 14, 
17, 28, 44, 58, 72, 84, 97, and 115. The settings of the 
system parameters are presented in Table 11^ The power price 
information is acquired from the PJM | |25[ data, and the 
aggregators are assigned with the prices of different buses. 

In the test system, 30 000 EVs are accommodated. We 
consider two models of vehicles, namely the Tesla Model 
S AWD-85D (Telsa Model S) and Nissan Leaf 2014 model 
(Nissan Leaf). We assume 60% of all EVs in the system are 
Tesla Model S EVs with 85 kWh batteries and 34 kWh/100 
miles energy expenditure performance, and the remaining 
40% are Nissan Leaf EVs with 24 kWh batteries and 31.64 
kWh/100 miles energy expenditure performance. The maxi¬ 
mum charging rates for these two models are 22 kW and 6.6 
kW, respectively, and the maximum discharging rates are also 
set to the same values | |26) . Among all EVs in the system, 80% 
EVs enabled bi-directional V2G operations. 5% of the EVs 
will perform a late departure up to a maximum of one hour. Eor 
practical situations, these parameters can be manipulated by 
the EV owners. The European Commission Strategic Energy 
Technologies Information System reported a mobility survey 
on the driving and parking patterns of European car drivers 
0- The results of the survey are utilized to formulate the EV 
driving dynamics. Similar methods have also been adopted in 
the literature, e.g, 0- 

We perform simulations on a horizon of 72 hours, and the 
length of each time slot is 15 minutes. Thus the total profit of 
all 288 time slots are combined as the performance metric. As 
Problem Q cannot be solved directly in a timely manner, we 
first investigate the performance of the distributed heurisric- 
enabled aggregator profit maximization strategy proposed in 
Section IV, whose optimality will be studied later. The dis¬ 
tributed solution is labeled “All” for ease of demonstration. 
In addition, four other variants are considered to demonstrate 
the relationship of Eq. 0 and LMP calculation, and their 
contributions to the profit maximization performance. These 
four variants are developed by removing one or multiple 
components from “All”: 


TABLE II 

Parameter Settings 


Param 

Value 

Param 

Value 

Param 

Value 

ginax 

23 

AR'' 

$0.015/6hrs 

„ch „dch 

0.9 

At 

15 minutes 

AR° 

$0.015/6hrs 

cmin 

0.0 

R'^ 

$0.08/kWh 


6 hours 

^max 

_ 

1.0 


$0.I0/kWh 


6 hours 

0.9 


TABLE III 

Profit Maximization Pereormance 


Mode 

All 

NoTrade 

NoLMP 

Planning 

Greedy 

Profit($) 

134295.1 

110269.8 

118534.4 

109750.6 

75678.9 


1) The NoTmde mode removes all possible trades between 
aggregators. This mode works like a V2G enabled 
version of | [22l . 

2) The NoLMP mode removes the LMP adaptation step. 
Energy trades among aggregators are allowed. 

3) The Planning mode removes both trading steps and the 
LMP adaptation step. This mode works similar to a 
V2G-enabled version of GD. 

4) The Greedy mode performs the greedy charging strategy. 
Upon the arrival of an EV, it is charged at the maximum 
rate until the SoC requirement is met. 

Table [represents the total profits generated by the compared 
charging approaches. It can be easily seen that our proposed 
approach can significantly increase the aggregator profits, 
and both the aggregator trading and LMP adaptation have a 
positive influence on the profit maximization process when 
compared with the Greedy mode. In addition, the computa¬ 
tional time of the proposed approach is also crucial as the 
algorithm is supposed to be implemented online. With 30 000 
EVs in the system and 24 optimized time slots, each iteration 
of our proposed algorithm can be finished in 9.13 seconds. 
All time slots can be optimized within a maximum of six 
iterations, i.e., one time slot can be finished in one minute. 

In addition to the profit maximization performance com¬ 
parison, we also investigated the optimality of our proposed 
approach in Section IV (“All” mode) comparing with the 
original one in Section III, and the behavior of aggregators in 
response to the power price. The detailed results are presented 
in f^ , and we observe that “All” is only 3.7% worse than 
the optimal solution of the convex-relaxed problem, and the 
gap between the true optimal and our approximation must be 
smaller. Eig. [^demonstrates the total charging power dynamics 
of our proposed approach with the changes on the average 
power purchasing price. We can observe that the charging 
power is mostly high when the price is relatively low. This 
shows the efficacy of the multiple time slot optimization in 
saving aggregators’ power cost. The impact of introducing 
the proposed strategy on the power system stability is also 
illustrated in | [2^ , and the conclusion can be drawn that 
adjacent buses to aggregators will be influenced more by the 
power line congestion, which is represented in the form of 
LMP changes. 

VI. Conclusion 

In this paper we propose a coordinated charging problem 
considering lET for maximizing the profit of multiple aggre- 












Fig. 3. Total charging loads for all time slots. 


gators. A model predictive control based centralized iterative 
approach is devised to find the optimal EV charging strategies 
for profit maximization. Considering the non-convexity nature 
of the optimization problem, we develop an aggregator-level 
coordinated charging heuristic to construct the EV charging 
schedules. To exploit the potential of employing energy trade 
among aggregators for profit maximization, we propose an 
auction-based heuristic to handle the trading details. In addi¬ 
tion, the iterative approach utilized by our developed strategy 
can further adapt the EV charging schedules to the grid 
congestion cost. To validate the performance of the proposed 
approach, we employ a ten-aggregator system with 30 000 
EVs for simulation. The system is installed on an IEEE 118 
bus system, and the simulation is performed on a 72 hours 
time span. The charging schedule developed by our proposed 
approach can create more profit for the aggregators than 
the compared strategies. Moreover, the proposed distributed 
optimization-based heuristic is compared with the relaxed con¬ 
vex optimization problem. The result shows that our heuristic 
can achieve near optimal performance. 
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